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Since the development of causal path analysis by 
Wright, both the standardlzea and unstandardi2ed prccedures in path 
coefficient estimates have been advocated in the related literature^ 
but the superiority of Dne technique over the other has not been 
demonstrated empirically or mathematically. In this study^ the 
sampling properties of both coefficient estimates in a chain 
prediction model are compared under known conditions by Monte Carlo 
methods at systematically selected parameter points which include 
sample rize and Var (E^sub-one, E--sub^two) * The implications^ 
interpretation^ and problems of both methods are discussed with 
reference to the evaluation of learning hierarchies, (Author) 
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A. The Backgr ound 

Since the developnient of path analysis by Wright (1918) , 
the technique of path analysis has gradually become popular in the 
medical I biological* social and political sciences* 

In the related literatures two types of path coefficient 
estlniates have been advocated, Tukey (1954) ^ Turner and Stavens (1959) 
have favored the use of unstandardized path coefficients, Wright (I960) 
restated some of the unique advantages of unstandardized procedure ^ but 
in general he favored the standardized form partly because of greater 
convenience in analysis. Applications of both types of approaches to path 
analysis have been reported. 

Although two methods of calculating path estimates are 
aviilablet the suptriorlty of one method over the other has neither been 
demonstrated mathematically nor empirically, Zellner (1971) reported a 
Mnnte Carlo study of the sampling properties of the two-stage least 
square beta weight in a tvyo-equation simultaneous equation model , but 
the only parameter he used was sample size and his empirical 
estimates were based on inforfnation from only 50 replications. Extending 
Zellner*s work, the overall purp^jse of this study was to compare 
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emp1r1cany under known conditions the sampling properties of both 
the unstandardized (t,e. derived from the beta weights) and the 
standardized (i.e. derived from the correlation coefficients) path 
coefficient estimates at selected parameter points. 

B. The Proble m 

Turner and Stevens (1959) identified three basic models 
of open path networks i namely: 

(1) the simple multiple prediction model, 

(2) the simultaneous prediction model » and 

(3) the chain prediction model. 

An example of the first model is shown in Figure 1. 




Figure 1 

In this and the fonowing path diagrams, the X's are the observed 
error free exogenous (independent) variables, and the Y's are the 
observed endogenous .(dependent) random variables. The arrows indicate 
the direction of influance. In intirpritlng thi path diagram, the 
variable at the head of one or more arrows is regarded as being a function 
of just those variables at the tails of these same arrows. The notation 

O 
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'a' indicates that the path coefficient is derived from the regression 
or partial regression coefficients, When the variables are standard- 
ized or when the path coefficient is derived from ths-correUtion or 
partial correlation coefficients » the notation for the path c«*fP»&tent 
in this paper 1s 'tt' . 

• ■ The model equation in Figure 1 is: 

(1) Yj = «^ + Xa + Xb + ■ , 

This is the familiar tnultlple linear pradfctlon model with two 
Independent variables and the path coefficients are the ordinary 

S 

partial regression coefficlants. In symbols T 

(2) ^ bi 

This model Is straightforward since the distribution of b^. (1n matrix 
notation). is 

(3) b - N { B , g2 (X'X)^^ } ■ 
where is the variance of e^. 

An eKample of the second path model, the siinultaneous prediction 
models IS shown in Figure 2. 




Figure 2 



correlated errors 
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The model equations in this case are 

(4) Yi - «^ + a^^ h + £i 

This model is also straightforward, since the distribution of b Is 

(5) b i H { B , I ( riv^ } 

where |_ in this case is a two-by-two error variance-covariance 
matrix. In the path technique as described by Turner and Stevens (1959, 
pp. 243 - 244)i the errops of the two dependent variablec.- (1 .ei- e^- and 
in ) were assumed to be uncorrelated. 

The third model. I.e. the chain prediction model 1s more 
unusual. The path diagram of the simplest of this type is inustrated 
in Figure 3. 



h > h—^h 



\ '2 



Figure 3 

The model equations for this example are 

and the reduced regression equations are then found: 
(7) Yl = «^ + Xa + 



correlated errors 
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The solution of the above path estimates can be approached by non- 
linear least square techniquess but due to the computations involved^ 
linear least square technique is generally used* In terms of ordinary 
regression coefficients, the path estimates can be found as: 
(8) • a^^ ^ 



«21 



b 



la 



where a,, is distributed as 



o2 



(9) N { . " } 

. SSx . 

but no comfc-'table distribution theory is ivanable for iji since It 

is a ratio of ^two correlated random variables* Should both a , a b 

la 21 

expressed in the standardized form* i.e, 

(IC) - 
and - s ^2a 

^la 

then naither tt^ nor jr^^ would have a comfoftable distribution 
theory. It is partly because of this reason that Tukey (1954), and 
Turner and Stevens (1959) adyocated the used of the unstandardizid 
path coefficient fstirnatas.- Our Qonjecturei were that both and 

ffji would have standard errors bigger than that' of S and ? . 

' ' la M' 

respectively and also that there might be some other undesirable 
sampling properties. Thus the specific objective' of this study was^ 
to to Investigate, by the Monte Carlo method, some of the empirical 
sampling behaviors of a^i and ynder selected parameter points . 
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including sample size and differences in the error \/ar1ance-covariance 
matrix of the dependent variables - 

C, Methodology 

Before the collection of empirical data, some of the 
theoretical values had first to be solved. In the following notations * 
^ is an error free variable* a and it are the population parameters. 
According to the model (see Figure 3), 
(11) ^ 



^Xa 



Y2 



21 



By specifying the parameters of , aia\ aa, 0^21,^2 and gC^^e^), 
the theoretical varecovarlances of X^* Yj and Yg can be derived to be 



(12) 




X 

a 


Xa 


^1 




h 





of .Cf2 + £j2 

1* Xa =1 

J\BL 



4l «fa 4z 



Based on this cgvarlance;. matrix, the population standardized path 
coefficients of ttj^ and Tr^j (1n its proper context of 0^ and' oji 
can be found as 

(13) ""la - 



'21 



'la 
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The procedures of compiling the data in the Appendix are briefly 
surrmarlzed in the following.',* 

(1) The vector of i,idapendent variable (X) with specified 
M and o was generated using Marsaglia's technique as described by 
Jansson (1966, pp. 182 - 185). It has been fOLind that this technique 
is reasonably fast and accurate. The piiudo-rdnyoni numbers of this 
subroutine have been tested for normality and randomness. The two 
vectors of correlated error terms with specified y's, o's, and p 
were generated according to the multivariate law as described in 
Oansson il%6, pp. 186 - 187). This method 1s general and can be 
extended to any numbir of dimensions. 

(2) The vectors of dependent variables (Yj» Y^) were 
computed as 

(14) h ^ °i + ^la' i ^ ei 

where Ej and are the vectors of error terms with correlation of p. 

(3) The sample variance-covarlance matrix frotn the generated 
Victors of X, Yl and Y2 was then computed and all the path coefflcent 
estimates of Figure 3 (i.e. «ja. Sgj' ^la' Hi^ calculated. In the 
expressions . . 

(15) = Ji_ 

"u ■ . 

and 



* Progranimed 1n Fortran IV by Peter T.K, Tam 
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it can be seen that as and r, approaches zero, the values of 

la la 

0L22 iT^^ can approach infinity. To avoid such unusual poss1bnities 
a sample of ,X, Yl and Y2 was accepted only if its computed | b^^ | 

was greater than BC and the corresponding | r^^ | greater than CC 
where 

(16) BC ^ (O.OS) ( cy|i / SSx )^ 

and 

CC - ( BC ) ( g| / ffj )^ 

\ " 

An example of the b, -s and r ^s that were accepted and rejected 

la la . 

was shown In Table 5 in the Appendix. 

(4). Steps (1) to (3) were replicated 200 times for each of 
the 15 different cycles of parameter values (see Tables 1 - 4 in 
Appendix) within each of the four different ievels of iampla sizes. 

The main computer program and Its subroutines have been 
carefully checked before Its execution on the CDC 6500 at the Florida 
State University Computer Center, The CPU time is approximately 
100 minutes and the CM 1s 50 K. . 

Results 

As mentioned previously, the reason for using the Monte 
Carlo approach is that no existing comfortable distribution theory 
is available for ql^i and ir^^. It seems that our conjectures were 
confirmed. All the summary data were tabulated in Tables 1 - 4 in 
the Appendix* Some of the interesting findings were graphically 
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illustrated in Figures 4-11 and discussed in the following: 
(1) Figures 4 and 5 indicated some of the effects of 
Sample size and error variance on the standard error of the 
path estimates. These graphs depicted the information of cycles 
i, 6, and 12 of the Appendix Tables. In Figure 4, the broken line 
indicated the amount of standard error of Sg^ across the four levels 
of sample sizes. The empirical and theoretical standard error of 
aia are represented by the continuous and the red line respectively. 
It can be seen -that the empirical and theoretical values are almost 
identical. Figure 4 is divided in three parts, namely A, and C, 
each of which has error variance ( cj ) of 50, 400, and 900 .respectively. 
As error variance ( ) increased, the standard error of the estimates 
increased. But within each level of error variance, the standard 
error decreased with the increase of sample size. However, irrespective 
of the sample size and values of ei, the standard error of Sgi is 
larger than that of This fact generally holds under other 
parameter points as indicated in the Tables in the Appendix. 

(2) As illustrated in Figure 5, the sama conclusions of 
above are generally applicable to that of ^^j. In terms of sta. dard 
error, it's indecisive of the relative superiority of and t^^. But 
both have empirical standard errors larger than that of fi^^ and ff^^. 
It seems that except when the sample size is big or the. error variance 
small , the standard error of S^i or may be too big to be tolerable. 

(3) Figures 6 and 7 illustrated the effects of sample size and 

ei on the empirical skewness of the path estimates. However, such effects, 
if there is any, are not obvious in this study. Figure 6 compared the 
skewness of with that of a^, and in this respect, two points might 
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be fairly obvious. Firsts the skewness of agi was in general larger 

than that of a^^. Second ^ the skewness of was In general large 
and positives although in some other cycles as indicated in the Tables 
in the Appendix, the skewness of 3^^^ was large- and negative. As 
illustrated in Figure 7, the same conclusions were true with respect 
to the skewness of ft^^. Again the superiority of or tt^^ was 
indecisive but neither one of them was desirable in terms of the 
skewness of the estimates. 

(4) Figure 8 and 9 compared S^^ and tr^^ in terms of the 
empirical bias of the estimates. The amount of bias in general 
decreased when the sample size was large or when was small or 
both* Again the superiority of either 5^^^ or was not obvious 

' A A 

but both were inferior when compared with that of ^nd 
respectively* 

(5) One of the assumptions made by Turner and Stevens (1959* 
p*241) in the solution of the path estimates was that the covarlance 
among the error terms ( ) was zero. Although in our experiment 

we had Introduced various levels of error covarlance (1,e* aCeie^)). 
results of this present study indicated that, under the parameters we 
Investigated^ there were no definitive trends between error covariances 

A A 

and the sampling propertiis of otji and irgi . This fact was rnustrated 
In Figures 10 - 12. 

In conclusion to the above, it seems that the results of 

A 

this study do support our conjecture that both a _^ and ^ in the 
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chain prediction model (see Figure 3) have indeed some Linusual 
sampling properties /Sreater cart in the 1ntirpretat1on of the 
coefficients In path analysis should be exercised when the chain 
prediction model is involved. This point is particularly signlfiCQnt 
s:nce most of the reported path networks in the related literature 
do have this chain prediction component. 
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Figuri 4 , Efftcts of sample siii on the standard error of the unstdndirdUed 
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Figure fl , Effects of sample size on bias of unstandardl zed path estimate.s 
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Figufi 1 . Effects of sample size un bf^s of stfln<lard1zid pjlh estimates 
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Figure 12. Effects of error cevariance on the bias of (A) 
and s^i (0). information based on cyciii 5 - S. 
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Sunanary of statistics for sample size of 1,000, 



25 



o o in 

Q^ O CO 



o% O eN 



00 



1^ 



=i Q% O 



in o ^ 



o in ^ 



o in 
in ^ 



in cn 



in m 



eg 
- u 

^ CM #^ 

t? o o 



tn 



cn 



in p 
cn cTk 



CM 



WD 



m 



cn o 



^ C4 

N N M 



■M 




s 








s 

• 


• 




in 
<s 


m 


es 


S% 




en 

o 


rn 

CM 




CN 

en 


cn 


fM 



in 

c> o o 



CM 



^ a^ ?M cn 
o o o o 



cn in «n tm 

•=H t=l O CM 



m cn o 

i=H O CM 



en 

i=i ^© cn 



rn cn 



yi ^ fM 

o o o 



in 



«^ 3-4 



in 



CM 



■-i CM rH CM 











€h 
j' 


?M 


€n 

r 




1 




In 
1 


£ 
« 


CM 


O 


cn 
1 


m 


CM 


m 

CM 


00 


cn 


OS 


O 


OS 
CM 




r 


1 







cn i-H 



wn 



CM 


cn 


o 




S 


o 


o 


o 


l' 








CM 


•<t 








o 
















rn 


CD 




o 


o 




C 











cn 








CM 




cn 





I I 



*M cn 



cn 
O 



CM 




CM 


GO 


CM 


fM 




CM 




m 




• 



*n 








cn 


o 


















cn 




r—i 


€M 


o 


cn 








• 


CD 


CM 


© 






CM 


CM 


cn 






• 


4 





GO 


€M 








CM 


CM 












r 


t' 






CM 








o 










• 





eM 

O 



r 



O 



i 



1 1 1 




in 












lA 










1 




l" 










•-H 


1—1 






<M 






< t= 


< t= 











o 



# 

I 



fM 



en 
o 



r 



CM 

o 



c5 



o 



o 



r 



o 



o 



€3 



< C3 



p 



i 

■ ^ 



«^ *H Iif5 j 



! 



' ^ \s\ m (- - - ■ 



I *N T»i ^, ry 



H N I 



ff» ^ e 
w «o 

^ PJ ff! 

«e ^ 



II ^ rr 

U iSJ U> ff! 



m ^us <^ 4^ r4fn as ^ m ^ 



fu *M *^ 



^ nQ m & 
» *=i ffk 



iSi *H QTi nQ 
^ CT* m ijp fii 

fVJ IM «^ ^ I 



N, >4 



J ^. , 

H ^ <r* i 



s aj 
t5 tfi K ^ K 

^ ^ fy ^ 5^ 
fn \r% tn i£3 ^ 



(t^ @ CT ^ N. 



^ yi 

= ^ ^ . - *^ d' J _ , — . , w ^ 

^fUfSuiPn f ^ i£% tn ^ t-^ iii ^ .* 1^ 



_ _ _ , fg ^J 



! 

m -4 i^'m ^ ^ * liMsJ *iJ If* s J N ^ yi ly ,1. 



Id 
i 



*B ^ ^ €"1 €f* 4£ 

J fd ^ ^ ftj Ui r*j I 
iis ivi H 4-1 rx- f<>. ( 



1 I 

Ufi Pj ^ B P9 J| 5^ Ui S ^ ib 
iFl lii is^ l/i -ri 4j 



»| f^i ^ i-t 



CP (X* x£5 

p-j ,n 



«^ s r» 
• • ■ 



li 

in 

Vi 



Of . 



fvj 

!^ f.. .i 
-t ^ ^ 



ft VSi ^ 

N. ffi i£l ^ 

■ ■» * • v i 

.4 H ftp H 



ipk i\i 
'^1 pfi 



j\ a- 

fSJ rf 
r*j M M 



f rii J- 

1^ ly i*j f\i Pi u 

* • • • • L_ 

fu ftj H i M H 



N *u r\j 

^ t4 Pii TJ i 

S flf^ ¥t rTi 

!y O n K IV 

H *i N *■* f^i T^i 



N P5 j- ^ K Pi r-^ fn 



fti 04 ^ iiS j 

rj H i(> ^ 

^ i/) * rf 

^ <9 

T( ^( T« pi j 



rp ff, 'f) ,^ 

•j< ifv m r*> K 

T »=# ^ ip f,s, 

H F< {\J H 



*^ ♦ -i f 



^" ^ Ui l*J N. 

J3 I fn K- ir» 
« H B 
r r s ' 

'9*99 

Pi *H pj pj 



^ ^ 
iP H 

^ f»j if* 

^ m 



UJ 


«M i 






CL 




Ui 






N= di y ^ 




^ in Pj 




hs. (\i 




P3 PJ 




■ ' • ■ 















fy oi fSi B ^ 
a a\ 

iN fy .4 j= ^ 



H Pj 

, s Pi p= , 

*^ B m ^ LP 



hs. ^ ii^ ^ ^ ^ ^, ^ 

1^ ^ ifs yj s r^t >[j 
^* ^ !P d fP ip ^ (O 



§ W J- 

Pi ^'i f Mi ffl 

f^| ifi T=| 

fv .4 u« I'J 

fy H iJ ~ H 



E» pu r** 



PJ P,i ^ ^ 



^ 0^ ^7 Sr « 

N i» pJ *e ff* 

J N. fVJ H « 

ry ^ ^ IP ^ 

* * • e • 



N & ^ ^ 

^ ^ hs # 
tn ijj \s\ 
*N tw J' 

• p ■ • « 

Pj PJ i W4 ^ 



m €\i ^ 

us \£% ^ 

^ pj li^ ^ M 

<<i q& > «^ K N 
• • • • > i 

^ pi 9^ ^ cv 



D e« ft ^ 
S ^ 1-1 ^ 
s# ^ tn m in 
t ^ l-i yj 

fSi 



yj ^ ^ 

^ t4 iA 

@ H ^ ^ 

PJ CM iS pJ pi, 
■ • ■ ■ 
t4 



Pi Ps. P3 
& iP ^ 

S ^ PMP 



{ ' 

!'* tf\ (fl in If H rl £ ? 

N, fp m ifi Pi m j\i ^ ir 

'*( (P !"i *4 F ^ f*^ pj f" 



S ^ 1! 0 * .n 4ft *S .#1 
f^' fft Pi ^ ^ in PJ p^ 

IH ? ^ Lf. NJ ^ P?J J ; 



-n r? ft 

£ * # *=l J 
J 



■H pj . 

ise fn 
PJ Pi ^ ! 



'p ift IP 

f J s3 tH Ml 



Pj N, 4 

M P 

ick ^p ' 



a* 

p^ p* N ,. 
*e IP 



*P JP {} p<a 

PP Pj .f 4?* .f *M €1 PJ ^ F 



tvj 'p I ^* Pi 



JJ -ij J , 



S (ti ^ 
© 



UJ 



is 



I 



N. LP pJ 

to %fi 1^ 
N» H tf( 

M -T e 



-71 H iS IP ff 



ifi ^5* isi 

10 T-i f ^ 

H ^ # 
iii at ^t' <i 



i t ^ N. ?n B 1^ pn ^ 
PJ Pj in fn 'P f a p^ rfi 

l'^ «P !f 1= PJ P. **' nj 



t4 CVI P5 

Ift f\j ^ 

€j 4> fSi iSf f»7 

^ IN PJ j= 



m m ui ij 
«H us ^ ^ H ^ 
^ ^ ff» H PJ 1^ 
PJ ^ sH ir 



H IP PJ s P5 

•/^ w s <s 

0^ j ^ ^ T- 

y3 ?H !^ fU ^ 



^ ^ © 

»S pa 

«^ 4# 

(y H 



t*7 ^ r' 

«^ PJ lA pa ^ 

P- N5 fi, 

4© ^ PJ PI ^ 



^ i 
\D fn ^ t 

^ ^ ^ c 

t3 ^ 5 



to 

31^ 



M fn 

MS 



Of <3 



er|c 



